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ABSTRACT 

Generalized  time -reversed  space-time  block  codes  (GTR-STBC)  are  introduced  as  a  conceptual 
tool  to  examine  the  impact  of  unequal  power  allocation  in  aeronautical  telemetry  channels.  Two 
transmitting  antennas  are  employed  to  exploit  partial  channel  state  information.  GTR-STBC  are 
observed  to  perform  the  best  trade-off  between  the  signal-to-noise  ratio  and  inter  symbol  interfer¬ 
ence.  It  is  also  observed  that  the  optimum  transmitter  power  profile  for  the  measured  channel  is 
significantly  different  than  that  in  the  statistical  channel  model. 


INTRODUCTION 

Aeronautical  telemetry  deals  with  a  wideband  air-to-ground  communication  link  that  usually  com¬ 
prises  of  frequency  selective  fading  channels.  In  addition,  severe  size,  weight  and  power  limi¬ 
tations  are  imposed  on  the  airborne  transmitter.  Highly  power-efficient  transmitters,  such  as  ra¬ 
dio  frequency  (RF)  power  amplifiers  operating  in  full  saturation,  and  constant-envelope  modula¬ 
tion  schemes  are  employed  to  overcome  these  limitations.  Since  frequency  selective  fading  leads 
to  inter-symbol  interference  (ISI),  techniques  for  mitigating  ISI  using  the  aeronautical  telemetry 
standard  IRIG  106  [1]  (especially  the  shaped  offset  quaternary  phase  shift  keying  version  TG 
(SOQPSK-TG))  are  of  immediate  interest. 

The  goal  therefore  is  one  involving  single-carrier,  constant-envelope  modulation  operating  in  a 
frequency  selective  channel.  This  is  clearly  an  equalization  and/or  diversity  problem.  Diversity 
reception,  especially  by  widely  separated  antennas,  is  an  obvious  but  expensive  solution.  Another 
feasible  solution  is  transmit  diversity  when  the  size,  weight  and  power  constraints  allow  it.  The 
most  common  scenario  will  be  a  system  involving  multiple  transmit  antennas  with  a  single  receive 
antenna,  which  can  increase  the  reliability  or  throughput  in  multipath  fading  channels. 
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Under  flat  fading,  if  the  transmitter  knows  the  channels  between  each  transmit  and  receive  an¬ 
tenna,  then  spatio-temporal  coding  [3]  is  optimum  in  that  it  maximizes  signal-to-noise  ratio  [2]. 
If  the  transmitter  does  not  have  this  knowledge,  then  a  diversity-maximizing  orthogonal  design 
(such  as  the  Alamouti  code  [4])  is  optimum  [2].  Although  OFDM  is  generally  applied  on  a  per 
subcarrier  basis  in  frequency  selective  fading,  it  is  not  a  viable  approach  under  the  constraints  just 
described.  Approaches  suitable  to  single-carrier  modulation  are  of  interest,  specifically  the  time- 
reversed  space-time  block  codes  (TR-STBCs)  [5]— [9].  It  plays  the  role  that  the  Alamouti  code 
does  in  frequency  non-selective  fading  with  an  uninformed  transmitter.  If  the  transmitter  knows 
the  channels,  it  is  not  obvious  how  it  could  use  this  channel  state  information  since  SOQPSK-TG 
does  not  have  any  “variable  parameters”  such  as  the  number  and  location  of  constellation  points. 
The  only  other  variable  available  to  the  system  is  the  power  allocated  to  each  transmit  antenna. 

In  an  effort  to  explore  the  robustness  of  transmit  diversity  in  links  prone  to  multipath  fading,  chan¬ 
nel  sounding  experiments  involving  multiple  transmit  antennas  have  been  conducted  at  Caims 
Army  Airfield,  Ft.  Rucker,  Alabama,  [10,  11]  and  Edwards  Air  Force  Base,  California,  [12].  Initial 
experiments  involving  TR-STBC  applied  to  the  helicopter-to-ground  channel  sets  revealed  some 
curious  behavior,  which  was  described  in  [11]  for  QPSK.  This  same  behavior  with  TR-STBC, 
SOQPSK-TG  for  the  same  data  set  is  illustrated  in  Figure  1.  This  plot  compares  the  simulated  bit 
error  rate  (BER)  performance  of  TR-STBC  using  minimum  mean-squared  error  (MMSE)  equal¬ 
izers  to  the  simulated  BER  performance  of  a  link  using  only  one  of  the  two  available  channels. 
TR-STBC  performs  better  than  the  single-channel  link  using  only  channel  2,  but  worse  than  the 
single-channel  link  using  only  channel  1.  Clearly,  channel  1  is  better,  in  some  sense,  than  channel 
2.  In  fact,  channel  1  is  so  much  better  than  channel  2  that  incorporating  channel  2  into  a  TR-STBC 
system  only  makes  things  worse.  These  results  imply  that  there  are  cases  where  it  is  better  to  aban¬ 
don  traditional  TR-STBC  and  use  only  one  of  the  two  available  channels.  The  fact  that  this  curious 
behavior  can  occur  on  real  channels  prompts  one  to  ask,  “Can  this  curious  behavior  be  derived 
from  the  given  impulse  responses  of  two  channels?”  In  other  words,  it  appears  that  with  some 
form  of  channel  state  information,  it  is  possible  to  achieve  better  performance  than  TR-STBC. 

This  paper  answers  the  question.  As  a  conceptual  tool,  we  consider  a  fixed-power  transmitter  that 
allocates  a  portion  of  the  fixed  power  to  each  channel.  This  power  allocation  is  parameterized 
by  p,  the  proportion  of  total  power  allocated  to  channel  1.  Using  only  one  of  the  two  available 
channels  is  captured  by  the  case  p  =  1  (channel  1  only)  or  p  =  0  (channel  2  only).  Next,  we 
introduce  Generalized  TR-STBC  (GTR-STBC),  a  modified  version  of  TR-STBC  that  incorporates 
the  unequal  power  allocation  p.  We  also  develop  and  analyze  a  MMSE  equalizer  that  operates  on 
the  samples  of  SOQPSK-TG.  By  expressing  the  residual  mean-squared  error  at  the  equalizer  output 
as  a  function  of  p,  we  are  able  to  choose  p  to  minimize  the  residual  mean-squared  error  for  a  given 
pair  of  channel  impulse  responses.  This  criterion  neatly  captures  the  contributions  of  both  ISI  and 
additive  noise  at  the  equalizer  output  and  it  facilitates  the  derivation  of  mathematical  expressions 
at  the  equalizer  that  are  tractable.  Also,  generally  speaking,  reductions  in  mean-squared  error  lead 
to  reductions  in  bit  error. 

Our  numerical  results  show  that  the  mean-squared  error  criterion  is  capable  of  identifying  the 
cases  where  the  curious  behavior  occurs.  The  conceptual  tool  of  unequal  power  allocation  together 
with  GTR-STBC  define  a  simple  transmit  diversity  scheme  based  on  the  partial  knowledge  of  the 
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Figure  1:  Simulated  BER  plots  for  10-Mbit/s  SOQPSK-TG  using  MMSE  equalizers  over  a  pair 
of  representative  impulse  responses  measured  over  a  helicopter-to-ground  channel.  The  circular 
markers  are  the  BER  results  using  a  single  antenna,  the  square  markers  are  the  BER  results  using 
both  channels  with  TR-STBC. 


channel  by  the  transmitter.  Here,  the  transmitter  only  needs  to  know  p,  which  is  easy  to  compute 
at  the  receiver  and  send  back  to  the  transmitter.  This  simple  scheme  includes  transmit  selection 
diversity  (p  =  0  or  1)  and  traditional  TR-STBC  (p  =  1/2). 


MMSE  EQUALIZATION  OF  SOQPSK-TG 

The  system  considered  here  is  summarized  in  Figure  2.  As  before,  the  complex-valued  baseband 
equivalent  representation  [13]  is  used  for  all  signals.  Starting  with  the  block  diagram  of  Figure  2 
(a),  the  SOQPSK-TG  signal  xc(t)  is  transmitted  through  a  channel  impulse  response  hc(t)  whose 
output,  accompanied  by  thermal  noise,  forms  the  receive  signal  rc{t).  After  the  application  of  an 
anti-aliasing  low-pass  filter  with  impulse  response  ha(t),  T-spaced  samples  of  rjt)  are  produced 
by  an  A/D  converter.  Assuming  the  anti-aliasing  filter  does  not  distort  the  received  signal,  the 
samples  of  the  received  signal  may  be  expressed  as 


r{ri)  =  x(n)  *  h(n)  +  w(n) 


L2 

^  h(k)x(n  —  k)  +  w(n) 

k=-L  1 


(1) 


where 

r(n)  =  rc(nT)  x(n)  =  xc(nT )  h(n)  =  hc(t )  *  ha(t) \t=nT 

and  where  the  w(n)  is  the  n-th  sample  in  a  sequence  zero-mean  complex- valued  Gaussian  random 
variables  with  autocorrelation  function 

Rw(k)  =  ^E  | w(n)w*(n  -  fc)|  =  <J2w8{k).  (2) 
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v(n) 

(b) 


Figure  2:  A  block  diagram  of  the  system  that  applies  an  MMSE  equalizer  to  samples  of  a  received 
SOQPSK-TG  signal:  (a)  the  system  showing  the  continuous-time  signals,  the  anti-aliasing  filter, 
and  A/D  converter;  (b)  the  equivalent  discrete-time  system. 


Note  that  (1)  assumes  the  discrete-time  channel  has  support  on  —L\  <  n  <  L2.  As  a  first  step,  we 
apply  the  samples  r(n )  to  a  filter  matched  to  the  discrete-time  channel  to  produce  y(n): 

y(n)  =  r(n )  *  h*(—n) 

=  x{n)  *  h(n )  *  h*(—n )  +  w(n)  *  h*(—n ) 

" - - - v - '  V - v - ' 

/leq(n)  v(n) 

Leq 

=  ^  heq(k)x(n  -  k)  +  v(n) 

k=  Z/eq 

where  Leq  =  L\  +  L2  and  v{n)  is  a  complex  valued  Gaussian  random  sequence  with  zero  mean 
and  autocorrelation  function 

Rv(k)  =  | v(n)v*(n  -  &)}  =  ^Jhq{k).  (5) 

The  samples  y(n)  form  the  input  to  an  MMSE  equalizer.  The  MMSE  equalizer  is  an  FIR  filter 
with  coefficients  c(n)  for  —Lc  <  n  <  Lc  designed  to  minimize  the  mean  squared  error  between 
the  equalizer  filter  output  x(n)  and  the  sequence  x(n).  The  entire  system  may  be  represented  by 
an  equivalent  discrete-time  system  shown  in  Figure  2  (b). 

The  challenge  with  equalizing  samples  of  the  modulated  signal  is  that  the  underlying  continuous¬ 
time  waveform  is  not  wide-sense  stationary  [13].  This  fact  carries  over  to  the  samples  of  xc(t)  and 
is  captured  by  the  fact  that  the  autocorrelation  function  of  the  x(n)  is  of  the  form 

Rx(k,£)  =  {x(k)x*(£)}7  (6) 

that  is,  the  autocorrelation  function  is  a  function  of  both  sample  indexes,  not  the  difference  between 
them.  Consequently,  the  equalizer  filter  coefficients  are  a  function  of  the  sample  index  n.  It  is 
hard  to  see  how  this  solution  has  any  practical  utility,  especially  in  the  presence  of  a  real-time 
performance  requirement.  In  the  end,  the  designer  is  left  with  suboptimal  approaches  of  reduced 
computational  complexity  whose  accompanying  performance  penalty  is  acceptable. 


(3) 

(4) 
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The  simplest  suboptimal  approach  is  to  assume  the  signal  samples  are  wide-sense  stationary.  Here, 
the  autocorrelation  function  is  of  the  form 

Rx(k-£)  =  hi{x(k)x*(£)y  (7) 

that  is,  the  autocorrelation  function  depends  on  the  difference  of  the  sample  time  indexes.  The 
wide-sense  stationary  assumption  for  x(n)  greatly  simplifies  the  solution  as  the  equalizer  coeffi¬ 
cients  no  longer  depend  on  the  samples  index  n. 

The  question  is  now,  what  function  should  be  used  for  the  autocorrelation  function  Rx{k)l  Two 
approximations  are  investigated  in  [17].  The  first  is  an  empirically-derived  autocorrelation  function 
which  is  obtained  by  generating  a  large  number  of  samples  x(n)  and  using  the  standard  estimation 
technique  assuming  wide  sense  stationarity.  Given  L  samples  of  x{n)  for  n  —  0, 1, , , . ,  L  —  1,  this 
empirical  autocorrelation  function  is 


together  with 


Re(k) 


1 

2(L  -  k) 


L- 1 

x(n)x*(n  —  k ), 

n=k 


0  <  k  <  L  -  1 


Re(k)  =  R*(—k),  —L  <  k  <  0. 


(8) 

(9) 


The  second  approximation  is  to  assume  the  data  are  uncorrelated  which  generates  a  correlation 
function  of  the  form 


Ri(k)  =  <j2x8{k). 


(10) 


In  [17],  we  demonstrated  that  MMSE  equalizers  based  on  above  approximations  exhibit  almost 
identical  performance.  Consequently,  we  use  the  second  approximation  as  this  choice  simplifies 
the  computations  of  the  equalizer  filter  coefficients. 

Because  the  optimum  equalizer  coefficients  no  longer  depend  on  the  samples  index  n,  the  relation¬ 
ship  between  x(n)  and  the  equalizer  output  x(n)  is 


x(n)  =  c(n)  *  y(n)  =  c(k)y(n  —  k ). 

k=-Lc 

The  vector  of  filter  coefficients  that  minimizes  the  mean  squared  error 

1 2 


£  =  E 


x(n)  —  x\n) 


(ID 


(12) 


is  given  by 


c  = 


GGf 


Eb/N, 


H 


eq 


-1 


gt, 


(13) 
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where  c  is  the  (2 Lc  +  1)  x  1  vector  of  filter  coefficients,  G  is  the  (2 Lc  +  1)  x  (2 Lc  +  1  +  2 Leq) 
matrix 

^/eq(-^/eq)  *  *  *  ^eq(  -^eq) 

^'eq(-^eq)  *  *  *  ^eq(  -^eq) 

^eql-^eq)  *  *  *  ^eq(  -^eq) 

Heq  is  the  (2 Lc  +  1)  x  (2 Lc  +  1)  matrix  given  by 

he  q(0)  •  •  •  heq(—2Lc) 

Heq  =  ;  i  (15) 

heq(2Lc)  ■  ■  ■  heq(  0) 

where  it  is  understood  that  heq(k)  —  0  for  \k\  >  Leq. 

Vector  g  is  the  1  x  (2 Lc  +  1)  vector  given  by 

g  =  [Olx(Lc-Leq)  heq  0ix(Lc_Leq)]  ,  (16) 

where  Oix(£c_l  )  is  a  row  vector  comprising  Lc  —  Leq  zeros  (here  we  assume  Lc  >  Leq,  i.e.,  the 
equalizer  is  longer  than  the  channel),  and  heq  is  the  1  x  (2Leq  +  1)  vector  given  by 

heq  [(^eq(-^eq)  ’  "  ‘  he q(0)  '  '  '  he q(  Z/eq)j  .  (17) 

The  corresponding  mean  squared  error  is 

GG'  +  ^^  g'Y  <18> 


GENERALIZED  TR-STBC  (GTR-STBC):  NON-EQUAL  POWER  ALLOCATIONS 

USING  TR-STBC 

An  abstraction  (to  the  sample  level)  for  a  2  x  1  TR-STBC  system  is  illustrated  in  Figure  3.  Here 
the  system  transmits  2 L  samples  of  an  SOQPSK-TG  signal  sampled  at  N  samples/bit  over  two 
transmit  antennas  to  one  receive  antenna.  The  equivalent  discrete-time  channel  between  transmit 
antenna  1  and  the  receive  antenna  is  represented  by  the  impulse  response  hi(n)  for  —  Li  i  < 
n  <  L  i  2  whereas  the  equivalent  discrete-time  channel  between  transmit  antenna  2  and  the  receive 
antenna  is  represented  by  the  impulse  response  h2{n)  for  —  L2,i  <  n  <  L2,2- 

The  TR-STBC  encoder  partitions  the  sample  sequence  x(0), . . .  ,x(2 L  —  1)  into  two  sequences 
Xi(n)  and  x2 (n)  as  shown  in  Figure  3.  The  lcngth-2/.  packet  is  transmitted  in  two  intervals1  each 

'in  a  practical  implementation,  a  guard  interval  at  least  as  long  as  the  longest  channel  impulse  response  must  be 
inserted  between  the  two  intervals.  Here,  such  an  interval  is  assumed,  although  we  won’t  complicate  the  notation  to 
make  this  explicit. 
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amplitude  gain 


Vi(n 


2/2  (n) 


Figure  3:  A  block  diagram  of  the  TR-STBC  system  based  on  unequal  power  allocation  using 

0  <  p  <  1. 


spanning  L  sample  intervals.  During  the  first  interval  :/•](()),. . . ,  x  \  ( L  —  1)  is  transmitted  from 
antenna  1  whereas  x2(0), . . . ,  x2(L  —  1)  is  transmitted  from  antenna  2.  During  the  second  interval, 
x2(L  —  1) , . . . ,  x2 (0)  is  transmitted  from  antenna  1  whereas  —  x\ (L  —  1) , . . . ,  —  x\ (0)  is  transmitted 
from  antenna  2. 


Power  division  using  0  <  p  <  1  is  accomplished  by  modified  the  TR-STBC  system  along  the  lines 
illustrated  in  Figure  3.  Amplitude  scaling  is  applied  to  the  signals  entering  each  channel  so  as  to 
divide  the  power  between  the  channels.  Here  p  represents  the  proportion  of  total  power  allocated 
to  transmit  antenna  1.  The  traditional  TR-STBC  system2  is  a  special  case  for  which  p  =  1/2.  The 
square -root  is  used  in  Figure  3  because  the  amplitudes  are  what  are  being  modified — the  energy 
(or  power)  is  the  square  of  the  amplitude. 


The  received  signal  r(n)  is  given  by 

r(n)  =  y/p  xi(n)  *  h\(ri)  +  yjl  -  p  x2(n)  *  h2(n)  +  w{n),  (19) 

where  w(n)  is  a  complex- valued  Gaussian  random  sequence  with  zero  mean  and  autocovariance 
function  given  by  (2).  The  TR-STBC  decoder  partitions  r(n)  into  ri(n)  and  r2(n)  as  follows: 


ri  (n)  =  r(n )  for  0  <  n  <  L  —  1 
r2(n  —  L)  =  r(n)  for  L  <  n  <  2L  —  1. 


(20) 


These  two  sequences  are  given  by 

ri(n)  =  y/p  xi (n)  *  hi(n)  +  yjl  -  p  x2(n)  *  h2(n)  +  Wi (n)  (21) 

r2(n)  =  yfp  x*2 (-n)  *  hi(n)  -  yjl  -  p  x\(—n)  *  h2(n)  +  w2(n),  (22) 

2In  the  traditional  TR-STBC  system,  p  =  1/2  is  included  in  neither  the  development  nor  the  notation  because 
the  same  power  is  assumed  to  be  applied  to  each  channel.  Hence  there  is  no  need  to  account  for  it,  other  than  in 
normalizing  the  noise  variance. 
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where  w\ (n)  and  w2 (n)  are  related  to  w(n)  the  same  way  ri(n)  and  r2(n)  are  related  to  r(n).  The 
TR-STBC  decoder  processes  ri(n)  and  r2(n)  using  a  bank  of  filters  based  on  the  channel  impulse 
responses  hi(n)  and  h2{n)  as  shown.  The  result  of  this  processing  is  a  pair  of  parallel  sequences 
yi(n)  and  y2(n)  which  may  be  expressed  as 


Ui(n)  =  r i (n)  *  yf p  h\(-n )  -  r2(-n )  *  y/l  ~  P  Mn) 

=  x\(n)  *  p  h\(n)  *  h*(—n )  +  (1  —  p)  h2(n)  *  h2(—n) 


htq(n) 


+  wi(n)  *  yfp  h\(-n )  +  w2(—n )  *  y/l  -  p  h2(n )  (23) 


vi  (n) 


and 


V2 (n)  =  r i (n)  *  y/l-  p  h2(-n)  +  r2(-n)  *  yfp  hi(n) 

=  x2(n)  *  (1  —  p)  h2{n)  *  h2{—n )  +  p  h\{—n )  *  hi(n) 

' - v - - 


heq{n) 

+  y/l  -  p  Wi(n)  *  h2(—n )  +  y/p  w2(—n )  *  hfn) .  (24) 

^ - V - ' 

V2  (n) 

These  equations  show  that  the  equivalent  composite  channel  for  non-equal  power  allocation  is 

heJn)  =  p  hi(n )  *  h\(—n)  +(1  -  p)  h2{n)  *  h2(-n ) .  (25) 

' - V - "  ' - V - ' 


m  (n) 


'/i 

Because  the  support  for  hfn)  is  —  L1;1  <n<  L12,  the  support  for  rp(n)  is  —  (L1;i  +  Li)2)  <n< 
(Lu  +  L1j2).  Similarly  because  the  support  for  /i2(n)  is  — L2j1  <n<  L2;2,  the  support  for  r/2( 
is  —  (L2ji  +  L2j2)  <  n  <  (L2j i  +  L2j2).  Consequently,  the  support  for  heq(n)  is  —  Leq  <  n  <  L 
where 


n 


Teq  —  max  i  +  Li;2,  L2ji  +  L2j2 

The  noise  sequences  tq(n)  and  u2(n)  are  complex- valued  Gaussian  random  sequences  each  with 
zero  mean  and  autocorrelation  and  cross  correlation  functions 

{vi(n)v*(n  -  A;)}  =  |u2(n)u*(n  -  A;)}  =  v2wheq(k). 


'eq 

(26) 


E  |ui(n)u2(n  —  A;)  |  =0. 


(27) 

(28) 


By  way  of  summary,  the  TR-STBC  system  presents  to  the  equalizers  the  sequences  y\(n)  and 
y2(n)  which  may  be  represented  by 


yi(n)  =  xi (n)  *  heq(n)  +  vi(n) 
y2{n)  =  x2(n)  *  heq(n)  +  v2(n) 


(29) 

(30) 
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where  heq(n )  is  given  by  (25).  The  noise  terms  v± (n)  and  v2(n)  are  uncorrelated  zero-mean  Gaus¬ 
sian  random  sequences  each  with  autocorrelation  function  (27).  A  pair  of  equalizers  operate  in 
parallel  on  yi(n)  and  y2 (n).  Because  the  the  noise  sequences  tq(n)  and  v2(n)  are  statistically 
equivalent  and  heq(n )  is  common  to  both,  the  pair  of  equalizers  operating  on  y,  (n)  and  y2 ( n)  are 
identical  as  long  as  ai(n)  and  a2{n)  are  statistically  equivalent  (the  usual  case).  Any  equalizer  can 
be  applied  here  (linear  or  non-linear,  with  or  without  noise  whitening)  with  the  usual  performance- 
complexity  tradeoffs.  Here,  we  apply  the  approximate  MMSE  equalizer  introduced  in  previous 
section  because  the  MMSE  equalizer  permit  a  mathematically  tractable  analysis  for  the  resulting 
mean-squared  error.  We  leverage  the  analytical  expression  to  find  the  value  of  p  that  minimizes  the 
mean  squared  error. 

The  equalizer  filters  of  Figure  3  are  identical,  and  the  vector  of  filter  coefficients  is  given  by  (13) 
and  the  corresponding  mean  squared  error  is  given  by  (18).  In  these  equations,  G,  g,  and  Heq  may 
be  expressed  in  terms  of  power  allocation  parameter  p.  Using  (25),  it  is  straightforward  to  show 
that 


G  =  pGi  +  (1  —  p)G2 

(31) 

g  =  Pgl  +  (1  -  P)g2 

(32) 

Heq  =  pHx  +  (1  —  p)H2 

(33) 

where  Gi,  gi,  and  Hi  are  formed  from  ry  (n)  the  same  way  G,  g,  and  Heq  are  formed  from  heq(n), 
respectively.  Similar  definitions  apply  to  G2,  g2,  and  H2  with  y2{n).  Making  the  substitutions  for 
G,  g,  and  Heq  gives 

c  =  IVT1  (pg!  +  (1  -  p)g2^  (34) 

and 

£min  =  (T2X  1-  (pgl  +  (1  -  p)g2)  M"1  (pgx  +  (1  -  p)g2)f  (35) 

where 

M  =  (pG,  +  (1  -  p)G2)  (pG,  +  (1  -  p)G2)f  +  AE-  (pH2  +  (1  -  p) H2)  .  (36) 

Equation  (35)  is  the  desired  relationship:  for  a  given  pair  of  channels  hi(n)  and  h2(n)  ,  it  expresses 
the  mean  squared  error  at  the  output  of  the  MMSE  equalizer  as  a  function  of  the  power  allocation 
p.  Thus,  for  a  fixed  pair  channels,  one  can  choose  the  power  allocation  to  minimize  the  achievable 
mean- squared  error. 


NUMERICAL  RESULTS 

The  forgoing  analysis  was  applied  to  a  helicopter-to-ground  radio  link  using  190,716  channel  im¬ 
pulse  responses  from  the  Black  Mountain  flight  path  at  Edwards  AFB,  California,  captured  during 
the  channel  sounding  experiment  described  in  [10].  For  the  modulation,  we  assume  10  Mbit/s 
SOQPSK-TG  with  the  equalizers  operating  at  2  samples/bit.  In  each  of  the  above  cases,  we  com¬ 
pute  the  optimum  power  allocation  parameter  p  to  a  pair  of  channels  labeled  hi(n)  and  h2(n).  The 
physical  meaning  of  these  channels  depends  on  the  scenario  as  described  below. 
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p  p 

Figure  4:  Optimum  power  allocations  in  the  mean-squared  error  sense  for  the  measured  channel 
impulse  responses  from  the  Black  Mountain  flight  path  at  Edwards  AFB,  California,  using  the 
natural  normalization:  (left)  Eb/N0  =  10  dB;  (right)  Eb/N0  =  20  dB. 


Two  normalizations  are  applied  to  the  channels:  the  natural  normalization  and  the  equal-energy 
normalization.  Let  /iljU(n)  and  h2u(n)  be  unnormalized  channel  impulse  responses  for  the  two 
equivalent  discrete-time  channels  obtained  directly  from  the  channel  sounding  data,  and  let 


L  1,2  L2,2 

Ei=  Y  IM™)I2  e2=  Y  1^2, uMP 

n=— Li.i  n=—L2:i 


(37) 


be  the  energies  in  two  channels.  The  natural  normalization  uses 

hi(n)  =  ~^=hitU(n)  h2(n)  =  ~^=h2:U(n) 


(38) 


where  E  =  m&x{Ei,  E2}.  This  normalizes  the  stronger  of  the  two  channels  to  unit  energy.3  We 
call  this  the  natural  normalization  because  in  real  multi-antenna  scenarios,  it  is  often  the  case  that 
one  of  the  channels  is  stronger  than  the  other. 

For  the  equal-energy  normalization,  we  use 


hi(n) 


h2(n ) 


(39) 


Here,  both  channels  are  normalized  to  unit  energy.  This  is  more  typical  of  statistical  or  mathemat¬ 
ical  models  of  multi-antenna  propagation. 

The  numerical  results  were  produced  as  follows.  For  each  of  the  190,716  pairs  of  channel  im¬ 
pulse  responses,  the  impulse  responses  were  normalized  using  one  of  the  two  procedures  described 
above.  The  value  of  p  that  minimizes  (35)  for  Lc  =  5  x  Leq  [18]  was  computed.  The  results  cor¬ 
responding  to  the  natural  normalization  are  shown  in  Figure  4.  Here,  the  ISI  is  such  that  the 
equalizer  is  able  to  reduce  the  residual  ISI  to  the  point  where  signal-to-noise  ratio  is  the  dominant 
factor.  The  end  result  is  that  the  ISI,  signal-to-noise  ratio,  and  power  allocation  tradeoff  becomes 
mostly  a  tradeoff  between  signal-to-noise  ratio  and  power  allocation.  This  tradeoff  is  such  that  the 
mean  squared  error  is  minimized  by  applying  100%  of  the  available  power  to  the  lower  antenna 

'The  motivation  for  unit  energy  is  for  scaling  the  noise  variance  to  define  the  signal-to-noise  ratio. 


10 


80 

8  60 

c 

<d 

§  40 

8 

o 

5s  20 
0 


□ _ 

— - n_. 

0.2 


0.4  0.6 

P 


0.8 


Figure  5:  Optimum  power  allocations  in  the  mean-squared  error  sense  for  the  measured  channel 
impulse  responses  from  the  Black  Mountain  flight  path  at  Edwards  AFB,  California,  using  the 
equal-energy  normalization:  (left)  Eb/N0  =  10  dB;  (right)  Eb/N0  =  20  dB. 


[corresponding  to  h2(n)]  in  the  vast  majority  of  the  cases  (152,246  or  80%  for  Eb/N0  =  10  dB, 
130,989  or  69%  for  Eb/N0  =  20  dB).  With  the  equal-energy  normalization,  the  signal-to-noise 
ratio  penalty  associated  with  hi(n)  is  removed  and  the  tradeoff  reduces  to  an  ISI  and  power  alloca¬ 
tion  tradeoff.  In  this  case,  summarized  in  Figure  5,  the  situation  reverses.  The  mean  squared  error 
is  minimized  by  applying  100%  of  the  available  power  to  the  upper  antenna  in  the  vast  majority  of 
the  cases  (147,738  or  77%  for  Eb/N0  =  10  dB,  155,786  or  82%  for  Eb/N0  =  20  dB). 

The  channel  between  the  upper  antenna  and  the  receiver,  h\{n),  is  characterized  by  lower  signal- 
to-noise  ratio  and  lower  ISI  than  h2{n),  the  channel  between  the  lower  antenna  and  the  receiver. 
Using  the  natural  normalization,  the  equalizer  is  able  to  remove  ISI  to  the  point  where  signal-to- 
noise  ratio  is  the  dominant  factor.  In  this  case  applying  100%  of  the  available  energy  to  the  lower 
antenna  is  the  optimum  approach  in  most  of  the  cases.  Using  the  equal-energy  normalization, 
differences  in  signal-to-noise  ratio  are  removed  so  that  ISI  becomes  the  dominant  factor.  Because 
the  ISI  on  the  channel  from  the  upper  antenna  is  so  much  less  than  the  ISI  from  the  lower  antenna, 
applying  100%  of  the  available  energy  to  the  upper  antenna  is  the  optimum. 

It  is  interesting  to  compare  these  results  with  what  might  be  inferred  from  using  a  simple  statistical 
channel  model.  To  do  so,  we  use  a  simple  Gaussian  model  for  each  channel  such  as  that  used 
in  [19].  In  this  experiment,  channel  1  comprises  3  independent  and  identically  distributed  (i.i.d.) 
zero-mean  complex-valued  Gaussian  random  variables  and  channel  2  comprises  6  i.i.d.  zero-mean 
complex- valued  Gaussian  random  variables.  These  numbers,  3  and  6,  are  the  average  lengths  of 
hi(n)  and  h2(n),  respectively,  in  our  measured  data  set.  The  channels  were  normalized  using  the 
equal-energy  normalization  described  above  and  190,716  independent  realizations  were  produced. 
The  results  are  summarized  by  the  histograms  in  Figure  6.  The  temptation  is  to  think  of  the 
optimum  p  as  a  normally  distributed  random  variable,  but  it  should  be  kept  in  mind  that  this  is  not 
the  case  because  0  <  p  <  1.  For  modest  values  of  Eb/N0  we  observe  that  the  mean  value  of  the 
optimum  p  is  about  0.5.  Given  the  fact  that  p  =  0.5  corresponds  to  traditional  TR-STBC,  we  see 
that  the  simple  statistical  model  suggests  that  traditional  TR-STBC  is  the  best  on  average.  This  is 
in  contrast  to  the  conclusion  drawn  from  the  measured  channel  data,  where  a  strong  preference  for 
transmit  selection  diversity  is  observed. 

As  Eb/N0  increases,  the  optimum  value  of  p  appear  to  increase.  This  is  explained  as  follows:  as 
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Figure  6:  Optimum  power  allocations  in  the  mean- squared  error  sense  for  pairs  of  equal-energy 
channels  where  the  channel  coefficients  are  i.i.d.  zero-mean  complex-valued  Gaussian  random 
variables:  (left)  Eb/N0  =  10  dB;  (right)  Eh/N{)  =  20  dB. 


Eb/N0  increases,  the  contribution  to  squared  error  from  additive  noise  is  reduced  and  ISI  starts  to 
become  the  dominant  contributor.  The  optimum  p  >  0.5  means  the  system  prefers  to  allocate  more 
energy  to  channel  1  than  channel  2.  This  makes  sense  because  channel  1  is  shorter,  and  this  tends 
to  contribute  less  residual  ISI  at  the  equalizer  output. 


CONCLUSIONS 

A  criterion  that  allows  one  to  predict  when  it  is  better  to  use  transmit  selection  diversity  (i.e.,  one 
transmit  antenna)  or  the  diversity  achievable  through  TR-STBC  (i.e.,  two  transmit  antennas)  was 
developed.  The  criterion  is  the  residual  mean-squared  error  at  the  output  of  an  MMSE  equalizer. 
The  residual  mean- squared  error  was  not  only  a  mathematically  tractable  quantity,  but  also  an 
excellent  predictor  of  the  curious  behavior  illustrated  in  Figure  1 . 

The  above  criterion  was  applied  to  a  set  of  measured  channel  impulse  responses  collected  from 
the  Black  Mountain  flight  path  at  Edwards  AFB,  California.  For  each  pair  of  channel  impulse 
responses,  the  value  of  p  that  minimized  the  residual  mean-squared  error  (35)  was  computed  and 
used  to  form  histograms  to  summarize  the  results.  Transmit  selection  diversity  (i.e.  applying  all  of 
the  available  power  to  the  stronger  channel)  is  optimum  when  two  channels  have  unequal  gains  in 
a  2-transmit,  1 -receive  antenna  system  operating  in  a  frequency  non-selective  fading  environment. 
In  contrast,  we  concluded  that  on  a  frequency  selective  fading  channel,  the  optimum  approach 
was  to  apply  power  to  produce  the  best  trade-off  between  SNR  and  ISI.  The  optimum  value  of 
p  associated  with  the  GTR-STBC  system  described  in  this  paper  identified  this  trade-off.  There 
were  some  channel  pairs  for  which  p  =  0  or  1  is  the  optimum  (transmit  diversity  case)  and  some 
channel  pairs  for  which  p  =  1/2  (traditional  TR-STBC).  But  there  were  many  channel  pairs  for 
which  neither  of  these  is  optimum. 

The  results  on  our  measured  channels  also  indicated  that  transmit  selection  diversity  was  more 
common  than  traditional  TR-STBC.  On  a  statistical  channel,  such  as  the  one  used  in  [19],  tra¬ 
ditional  TR-STBC  is  the  best  thing  to  do  on  average.  This  is  in  contrast  to  the  results  from  the 
measured  channels.  Consequently,  the  optimum  power  allocation  in  a  real  setting  is  not  predicted 
well  by  simple  statistical  channel  models. 
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